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Abstract 

o ■ 

^sj . In the 1970's, Balas [2J [J] introduced the concept of disjunctive pro- 

gramming, which is optimization over unions of polyhedra. One main re- 
[•"t | suit of his theory is that, given linear descriptions for each of the polyhedra 

to be taken in the union, one can easily derive an extended formulation of 

f — ' the convex hull of the union of these polyhedra. In this paper, we give a 

generalization of this result by extending the polyhedral structure of the 
variables coupling the polyhedra taken in the union. Using this general- 
ized concept, we derive polynomial size linear programming formulations 

(~J , (compact formulations) for a well-known spanning tree approximation of 

Steiner trees, for Gomory-Hu trees, and, as a consequence, of the minimum 
T-cut problem (but not for the associated T-cut polyhedron). Recently, 
C^ . Kaibel and Loos [S] introduced a more involved framework called polyhe- 

dral branching systems to derive extended formulations. The most parts 
of our model can be expressed in terms of their framework. The value of 
our model can be seen in the fact that it completes their framework by 
an interesting algorithmic aspect. 
> 

1 Introduction 



Let Q :— {(x, y) 6 1R P x R 9 | Ax + By > c} be a polyhedron. The projection of 
Q onto the x-space is the polyhedron 

Proj x (Q) := {x e W | 3y G R 9 : (x,y) e Q}. 



Conversely, Q is said to be an extension of the projected polyhedron P := 
Jv> , Proj a ,(<5), and the system Ax + By > c is an extended formulation for P. The 

5^ ■ extended formulation Ax + By > c is compact if q, the number of rows, and the 

input length of each entry of the inequality system is polynomial in p. 

In the 1970's, Balas [2JH] introduced the concept of disjunctive programming, 
which is optimization over unions of polyhedra. Below we restate a well-known 
result saying that, given linear descriptions for each of the polyhedra to be taken 
in the union, one can easily derive an extended formulation of the convex hull 
of the union of these polyhedra. Recently, Kaibel and Loos [5] introduced a 
powerful framework called polyhedral branching systems that generalizes Balas' 
result as well as the framework of Martin, Rardin, and Campbell [TT] to derive 
extended formulations from dynamic programming algorithms for combinatorial 
optimization problems. 

In this paper, we consider a generalization of the concept of disjunctive 
programming whose polyhedral aspects are covered by the framework of Kaibel 
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and Loos [5] . The motivation of our model lies in the algorithmic interpretation 
of disjunctive programming as a two- level approach to solve a linear optimization 
problem. Since the generalization is straightforward, while the framework of 
polyhedral branching systems is quite involved (at least, it would need some 
space to explain it), we pass on a description of that framework and directly 
start our considerations with disjunctive programming, its interpretation as two- 
level optimization model, and its consequences for extended formulations. 

Given a finite collection of polyhedra P z , i £ X, where X is a finite index set, 
a disjunctive program is a mathematical program of the form 

max w T x 

s.t. x£ |J P\ (!) 

iez 

By a well-known result of Balas [3j, given complete linear descriptions of 
each of the polyhedra P % to be taken in the union, one can describe the convex 
combination of the union of the polyhedra by an extended formulation. 

Theorem 1.1 (Balas [3]). Given polyhedra 

P l = {x G R n | A l x > b 1 } = convV' + cone if, i e X, 



the following system: 



•£** = <), 

- Xib 1 > 0, i € 1, 

A,>0, i€l 



(2) 



provides an extended formulation for the polyhedron 

P x := conv |J V i + cone |J R l . 

In particular, denoting by P the set of vectors (x, {x l , \i\iez) satisfying ([2]), 

(i) if x* is a vertex of Pi, then (x, {x l , Xi}i^i) is a vertex of P, with x = x* , 
(x k ,X k ) = (x*. 1) for some k £ X, and (x\\i) = (0, 0) for i G X \ {k}; 

(ii) if (x, {x l , Xi}i£i) is a vertex of P, then (x k , X k ) — (x, 1) for some k G X, 
(x l , Xi) — (0, 0) for i G X \ {k}, and x is a vertex of Px- O 

By Theorem 11.11 the disjunctive program (fTJ) can be solved by solving the 
linear program max{w T x \ (x, {x l , Xi]n=.x) satisfies ([5])}, provided we are given 
linear descriptions of each polyhedron P l as required. 

From an algorithmic viewpoint, to solve (TTJ), one usually would compute an 
optimal solution of each subproblem max.{w T x \x G P 1 }, and then one would 
choose the best (or a best) among them. This two-level approach is reflected 
in the extended formulation. For simplicity, let us assume that, in the first 
phase, for each subproblem m&x{w T x \ x G P 1 } an optimal solution x l exists. 
Let A G R x be the (variable) vector whose components are the A^. Then, one 
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defines a vector w by w; := w T x z , i £ I, and solves, in the second phase, the 
linear program max{w T A | A £ A} over the simplex 



{A e R 1 : ^ A, = 1, A, > Vi e 1}. 



Given linear programs max{w T x | x £ P 1 }, it is, however, not always in- 
tended to optimize over the union of these polyhedra, but sometimes the sub- 
problems are part of a more complex optimization problem. 

As an example, let us consider (the polyhedral version of) the minimum 
spanning tree problem over the metric closure of a weighted graph. This problem 
has relevance for the approximation of the Steiner tree problem, which will be 
discussed in more detail in Section [2 We denote the node and edge set of a 
graph G by V(G) and E(G), respectively. Given a graph G with a nonnegative 
edge weighting w : E{G) — > R + , the metric closure of (G, w) is the pair (K, w), 
where K is a complete graph on node set V(G) and w(e), for e = {s,t}, is 
defined to be the length of a shortest path connecting s and t in G w.r.t. to w 
if there is any such path, and otherwise w(e) := +oo. The aim is now to find 
a spanning tree T of K minimizing w(E(T)) :— Y, e &E(T) ^( e )- The approach 
obviously consists of a two-level model. In the first step, we solve a so-called all- 
pairs shortest path problem, and in the second step, a minimum spanning tree 
problem whose input data are given by the output data of the first problem. 
For each of the two problems, there are known several linear programming 
formulations and, among these, even compact formulations. The question now 
arises whether or not these formulations can be brought together to provide 
a linear programming formulation for the entire problem. The answer to this 
question is surprisingly quite easy and can be given using a modification of (O . 
Let I :— E(K) be given by the edge set of K, and let min{^ w(e)x e | A e x > b e } 
be a linear programming formulation of a shortest s,i-path problem in G for 
each edge e = {s, i] £ E(K) = I. Then, replacing A by a linear characterization 
LT := {A £ R^^ : CX > d} of the spanning tree polytope, we propose the 
following model: 



min 


Y, W eXe 
eeE(G) 




s.t. 


x — Yl x£ 

e£E(K) 


= 




A e x e - X e b e 


>o, 




CX>d. 



e £ E(K), 



In any optimal solution (x, {x e , X e }e£E(K)), A is then a convex combination of 
the characteristic vectors of minimum spanning trees w.r.t. w : E(K) — >■ R + . 

The aim of the remainder of this section is to prove the correctness of this 
approach in general. 

For any polyhedron P = {x £ R™ | Ax > b} and any a £ R+, let P(a) := 
{x £ R™ | Ax > ab}. This implies that 

P = conv{«i,v 2 ,...,«fe} + cone{ri,r 2 ,...,r m } 
■& P(a) = conv{avi,av 2 , ■ ■ ■ ,avk} + cone{ri,r 2 , . . . ,r m }. 

Theorem 1.2. Given pointed polyhedra P l = {x £ R™ | A % x > b 1 }, i £ I and a 
0/1-polytope II = {A £ R 1 1 CX > d}, that is, II = conv V for some V C {0, l} 1 . 
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(i) (x, {x l }i£i, A) is a vertex of the polyhedron Q defined as the set of vectors 
(x,{x l }i(zi, A) satisfying 

x- J2 x l = 0, 

AV - XiW > 0, iel, ( 3 ) 

CX>d 

if and only if X is a vertex of II, and for each i £ X, x % is a vertex of P % 
if Xi = 1 and x l = otherwise. 

(ii) For any w £ R", consider the linear programs 

max{u; T x | (x, {x l } ieI , X) £ Q}, (4) 

Wi :— ma,x{w T x \ x £ P 1 } (5i) 

for i £ X, and 

max{u) T A A £ II}. (6) 

Then, ^ is unbounded if and only if vdi = oo for some i £ X. Next, let ((J) 
be bounded. Then, (x, {x*}i g x, A) is an optimal solution of ^ if and only 
if X^ x 1 is an optimal solution of (l5l)) for each i £ I with \ > and X is 
an optimal solution of ((BJ- 

Proof, (i) To show the necessity, let (x,{x l }nzx,X) be a vertex of Q. First 
suppose, for the sake of contradiction, that A is not a vertex of II. Then, there 
are /i, v £ H and < a < 1 such that A = a/i + (1 — a)v. Let I 1 be the set of 
indices i £ T with Xi = 0. This implies that \Xi = z/j = for i £ T' . Define now 
vectors y % := z l := x % for i £ T' , y l := Qx 1 , z % := ~-x l for i £ X \ I', as well as 
y := Yliex V 1 an< ^ z := Siei z% ■ Then, one easily verifies that, on the one hand, 

a(y, {y l }iex, /«) + (1 - a)(z, {z l } ie i, v) = (x, {x l } ie x, A), 

and on the other hand, (y, {y l }i£X, n), (z, {z l } 2eI , v) £ Q. Thus, (x, {x l } ie x, A) 
is not a vertex, a contradiction. 

Next, suppose that (x, {x*}i 6 x, A) is a vertex of Q and A a vertex of II. The 
latter implies that A £ {0, l} 1 . Assume now that x 3 is not a vertex of P 3 
for some j with Xj = 1. Then, x 3 is the convex combination of two vectors 
y J , z 3 £ P 3 . This, in turn, implies that (x, {x l }i£i, A) is a convex combination 
of the two vectors obtained from (x, {x l }nzx, A) by replacing the vector x 3 by 
y 3 and z 3 . Moreover, assuming that x 3 ^ for some j with Xj =0, we see that 
x 3 is a ray of the cone {x £ IR™ | A 3 x > 0}. This immediately implies that, also 
in this case, (x, {x l }t e x, A) is not a vertex of Q, a contradiction. 

To show the sufficiency, suppose that (x, {x 1 }^!, A) is not a vertex of Q. 
Then, there are two different vectors (y, {y l }i^x, /i), (z, {z l }j S i, v) £ Q and 
< a < 1 such that 

Oi{y, {y l }iei, /") + (!- a )( z , {^}iez, ^) = (S, {x l } ie i, A). 

First, assume that \i ^ v. Since both vectors are in II, this immediately implies 
that A is a convex combination of (i and v, a contradiction. Consequently, we 
may assume that [i = v = X £ {0, l} 1 . However, since (y,{y l }i e x,tJ-) and 
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(z, {z l }i e xi v) have to be distinct, and since the vectors y and z are just sums 
of the vectors y 1 and z z , respectively, we conclude that y 3 ^ z 3 for some j G X. 
If Xj = 1, it follows that x 3 was not a vertex of P\ Finally, if Xj = 0, then 
ay 3 + (1 — a)z 3 — x 3 ^ or ay 3 + (1 — a)z 3 = 0, and hence P 3 is not a pointed 
polyhedron. In either case, this yields a contradiction. 

(ii) If (UJ) is unbounded, then there exists a ray (r, {r J }j e i, p) ^ of P with 
w T r > 0, which implies r ^ 0. Since II is a polytope, it follows that p = 0. 
Moreover, since r is the sum of the vectors r 7 , w T r* > for at least one i G X. 
Since, by definition, A l r l > 0, this implies that r l is a ray of P % and (J5T|) is 
unbounded. Conversely, if (|5il) is unbounded for some i G X, there exists a ray r 
of P l with w T f > 0. Define (r, {r 3 }j e z, p) by r := r* := f . r 3 := for j G I\{z}, 
and p := 0. Then, we conclude that (r, {r J } je x, p) is a ray of Q and u/ T r > 0. 
Therefore, (UJ) is unbounded. 

Next, suppose that ^ is bounded, which means that (JSTJ) is bounded for 
each i el. This, in turn, justifies the definition of w. 

For any i 6 X and any a G R+, x* is an optimal solution of max{w T x | x G 
P*} if and only if ax* is one ofmax{w T i | x G P l (a)}. Now let (5, {x J }j 6 x, A) be 
an optimal solution of (UJ). Assume that w T x* > w T (X~ l x l ) for some x* G P', 
with i G I and Ai > 0. Then, (x, {x 3 }j e x. A) G Q, where x l := A;x*, x J := x 3 
for j e I\ {z}, and x := Xiez 2 ' 1 - Moreover, w T x > w T x, a contradiction. 
Next, for any A G II and any optimal solutions x l of (J5TJ) for each i G X, we 
derive that (x, {x J }j e x, A) G Q, where x := X^ex 2 -* an ^ x% '■— ^^ f° r i £ X. 
Thus, w T x < w T x. Since 





w T x < w T x 


«=> 


YZi^W T X i <Ei6X^ T ^ 


<*=> 


Eiex™ Oi£') <E ie x^ T ^ 


<=> 


X; iGl AjW T x l < E^x^™ 7 ^' 


«=> 


Eiex^m, <Ej £ x^^ 


<^> 


u> T A < w T A, 



it follows that A is an optimal solution of ©. 

For the same reasons, if x % is an optimal solution of (J3TJ) for each iel and 
A is one of ([6]), then (x, {x J } jS x, A) is optimal for ([4]), where x := Eiex^*- '-' 



Theorem 1 1 . 2 1 can be generalized in several ways. Of course, it can be easily 
extended for the case that II is not a 0/1-polytope but any other polytope 
contained in R^. More important for the following applications is the case in 
which the polyhedra P l and II are given themselves by extended formulations. 
The following theorem only generalizes those results of Theorem 11.21 that are 
relevant for the following applications. 

Theorem 1.3. Given extensions 9 := {(A, p) G R 1 x R ? \ CX + Dp > d} of 
a 0/1-polytope II C R 1 as well as Q z = {(x, y) G R" x W* \ A l x + B l y > b 1 } 
of pointed polyhedra P l C R", i£l Moreover, for any w G R™, consider the 
linear programs 



max 


T 
W 








s.t. 


A 1 . 


v l + B 


V - Xib 1 
CX + Dp 


>0. 

> d. 



el, 



(7) 
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Wi :— max{w x \ (x,y) £ Q 1 } (8i) 

for i £ I, and 

max{u) T A|(A,^) £ 6}. (9) 

Then, J7]) is unbounded if and only if Wi = oo for some i £ X. Moreover, in 
case that Wi < oo for all i £ X, ({x 1 , y l }iex, A,/2) is an optimal solution of ([7]) 
if and only if (x l ,y l ) £ Q l , i £ I, are optimal solutions of (|8i|) and (A,/l) G O 
is an optimal solution o/ ([DJ. □ 

An application of Theorem 11.31 can be found in Pochet and Wolsey [13] for 
a variant of the classical lot sizing problem with constant production capacities 
in which the capacity in each period is an integer multiple of a basic capacity 
unit. 

In the following sections, we give two further applications of our model. 
Since both applications are based on compact formulations for spanning trees, 
we recall two well-known spanning tree formulations. 

Given a graph G, the spanning tree polytope PTrcc(G) is the convex hull of 
the characteristic vectors x E £ ~Rl EI > G ' of spanning trees T C G. Recall that 
the characteristic vector \ F °f an y edge set F C E(G) is a 0/1-vector with 
Xe = 1 if and only if e £ F. As it is well known, the spanning tree polytope is 
the set of all A £ IR^C 3 ) satisfying the nonnegativity constraints A e > 0, e £ E, 
the equation \{E(G)) = \V(G)\ — 1, and the inequalities 

X(E{G[S})) < \S\ - 1 for all S c V(G), 2 < \S\ < \V(G)\ - 1. 

Here, G[S] denotes the subgraph of G induced by S C V, and for any F C 
E(G), X(F) := J^eeF ^ e - This linear description of PTr e e(G) has an exponential 
number of inequalities. 

According to [3 [TU1 E)] , a compact formulation for PTr ee (G) can be given as 
follows. For each edge e £ E(G), we introduce a variable A e , and for each for 
each edge e £ E(G), each node u £ e, and each node v £ V(G), we introduce a 
variable \i e ,u,v In the extended formulation, the edge set of a spanning tree T 
of G will be represented by the vector (A T , /x T ) with A T = x Ei " T ^ an d fj,^ u v = 1 
if and only if e £ E(T) and v belongs to the same component of u in T — e. 
Then, PTrcc(G) is the projection of the polyhedron n defined as the set of all 
(A, jj) satisfying 






X(E) = \V\ - 


-i, 




t,v},u,v {^{u.v} ,V,U " 




for all {u, v} £ E, 


^e — 2-^d f^e,v,w " 


fo] 


- all e £ E, w £ V, 


v£e 






/ j l^{u,w},w 7 v ~ ■!• 




for all {u, v} £ E, 


GV\{u,v} 






A e >0 




for all e £ E, 


Me,f,iy ^_ v 


for all e 


£ E, v £ e, w £ V, 



(10) 



onto the A-space. Here, V := V(G) and E := E(G). 

For later reference, we need the following result which can be found in [ID] . 

Lemma 1.4. II is an integer polyhedron. Moreover, (A,/i) 6 II is a vertex ofH 
if and only if A is the characteristic vector of a spanning tree T of G' , and for 
each e £ E(T) and each u £ e, the vector {fJ- e ,u,v)v£V(G') * 5 the characteristic 
vector of the component ofT — e containing u, and for each e £ E(G') \ E{T), 
^e,u,v — for all u £ e, v £ V. □ 
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The following directed formulation is due to Wong [TS] and Maculan [S]. 
Let D be a directed graph, with node set V(D) and arc set A{D). For a fixed 
node r £ V{D), the r-arborescence polytope, that is, the convex hull of the 
characteristic vectors of r-arborescences, is the projection of the system 

v(A) = \V\-l, 

v a > v™ for all a £ A, w £ V \ {r}, 

a w (S+(r)) = l for all w £ V\{r}, (11) 

a w (S+(v)) -a w (d~(v)) = for &Rv,w £ V\{r},v^w, 

a™ > for all a £ A, w £ V \ {r} 

onto the i^-space, where V := V(D), A := A(D), and for any v £ V, 5 + (v) and 
S~ (v) are the set of arcs leaving and entering v, respectively. This formulation is 
motivated by the fact that an arborescence with root r contains an (r, w)-path 
for each node w ^ r. Let now G be a graph and D the digraph obtained from G 
by replacing each edge {u, v} of G by the arcs (u, v) and (v, u). Then, Pxree(G) 
is the projection of (fTTj) extended by the inequalities 

^{u,v} = v uv + v vu for all {u, v} £ E(G) (12) 

onto the A-space. 

The remainder of this paper is organized as follows. In Section [5J we pro- 
vide a compact linear programming formulation for a well-known spanning tree 
approximation of the minimum Steiner tree problem. In Section [31 we present 
a compact formulation for Gomory-Hu trees. Based on this formulation, we 
derive one for the minimum T-cut problem whose inequalities depend, however, 
on the objective function. Moreover, we discuss the relevance of these formu- 
lations for finding a compact formulation for the perfect matching polytope. 
In Section 01 we briefly summarize our findings and point out some interesting 
open questions. 

2 A compact formulation for the approximation 
of the Steiner tree problem 

Let G be an undirected graph and S C V(G). A Steiner tree for S in G is a tree 
T C G whose node set V(T) contains S. Given a cost function c : E(G) —> R+, 
in the Steiner tree problem for (G, c, S), one wants to find a Steiner tree TCG 
minimizing c(E(T)). 

Let (K, c) be the metric closure of (G, c) . By a well-known result of Gilbert 
and Pollak [B], if T is a minimum Steiner tree for 5, and M is a minimum 
spanning tree in K[S] w.r.t. c, then c(E(M)) < 2c(E(T)), where K[S] denotes 
the subgraph of K induced by S. 

For the all-pairs shortest path problem (for the computation of the metric 
closure) as well as for the minimum spanning tree problem there exist com- 
pact formulations. Hence, using Theorem 11.31 it is easy to derive a compact 
formulation approximating the Steiner problem. 

For any digraph D, we denote by V(D) and A(D) the node and arc set of 
D, respectively. Let G be a graph and c : E(G) —¥ R+. In what follows, let D 
be the digraph obtained from G by replacing each edge e = {v, w} £ E(G) by 
the arcs (v, w) and (w,v). Define c : A{D) — > R by c((v,w)) :— c({v,w}) for 
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(15) 



(v, w) E A(D). Then, a shortest s, t-path problem in (G, c) can be modeled as a 
minimum cost flow problem in (D, c). Consequently, for each edge e = {s,t} E 
E(K[S\), 

min Yl c(a)x e a 

aeA(D) 

s.t. x £ (S+(s))-x e (5-( S )) = l, (I3e) 

x e (5 + (v))-x e (5-(v)) = forallue V(L>)\{s,i}, 

< < < 1 for all a E A(D) 

is a compact formulation of the shortest s,i-path problem in (G,c). 

Combining the shortest path formulations (|13e[) for each edge e E E(K[S]) 
with a compact formulation for the spanning tree polytope P Tree (^[5']) defined 
on the complete graph on S, we obtain, by Theorem ll.3l a compact formulation 
approximating the Steincr tree problem defined for (G, c) . 

Theorem 2.1. Let G be a connected graph, let c : E(G) — > R+ be a cost 
function, and let S C V(G). Then, 

min ^ c(a)5^a;« (14) 

06,4(13) eG-E 

x e ((5+(s)) -x e ((5-(s)) - A e =0 for alle = {s,t} E E, 

x e {S+(v)) -x e {S-{v)) =0 /oraibey(D)\{e},ee£, 

<-A e <0 for alia eA(D), eE E, 

x e a > /or a/Z a e A(D), e £ E, 

(\,(t) satisfies (JTUJ) (or (\,v,a) satisfies ([IT]). ([HflO . 

where E := E(K[S]), is a linear programming formulation of size 0((\G\ + \S\ + 
(c) ) | <S"| 2 ) whose optimum is at most two times the optimum of the Steiner tree 
problem for (G,c,S), where \G\ := \V(G)\ + \E(G)\. 

Proof. Clearly, the above linear programs (fT4"|) subject to (ITU)) . (fT5|) and (fTT)) 
subject to (TTT1) . (|T^|) . (IT5|) are compact formulations of the spanning tree ap- 
proximation of the Steiner tree problem. So it remains to prove the correctness 
of the sizes of these linear programs. We only check the first one. Clearly, it 
consists of 0((\A(D) + \S\)\E\) variables and 0((\V(D)\ + \A(D)\ + \S\)\E\) 
inequalities. Since V(D) = V(G), \A{D)\ = 2\E(G)\, and \E\ = \E(K[S})\ = 
i|5|(|5|-l), the number of variables and inequalities is 0((\E(G) + \S\)\S\ 2 ) and 
C((|G| + IS'DIS'I 2 ), respectively. Moreover, each coefficient c({v,w}) is counted 
2\E\ times. Hence, the input length of the objective function is 0((c)\S\ 2 ). D 

3 A compact formulation for Gomory-Hu trees 

Gomory-Hu trees have many applications in a wide area of graph theory. To 
mention some examples, they provide compact representations of minimum s, t- 
cuts for all pairs of nodes s,t of a graph, they can used to determine minimum 
T-cuts, and they play an important role in the design of minimum cost commu- 
nication networks. 
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After some preliminaries, we first consider minimum-requirement trees that 
turns out to be the same as Gomory-Hu trees. As we will see, our approach 
directly leads to a compact formulation for Gomory-Hu trees. 

For any subset of nodes U C V(G) of a graph G, let S(U) be the set of 
edges of G connecting U and V \ U. A cut in G is an edge set of type S(U) 
for some ^ U C V(G). U and V(G) \ U are called the shores of £(J7). For 
two distinct nodes s and t of G, an s, t-cut is an edge set of type S(U) such 
that s € U and £ G V(G) \ U, or vice versa. We sometimes write 5g(U) instead 
of 6(U) to indicate that Sg(U) induces a cut of G Given a capacity function 
c : E(G) —> R+ and a cut if C E(G), the number c(-ftT) is called the capacity of 
if. The minimum s, t-cut problem asks for an s,£-cut K C E{G) of minimum 
capacity. 

As it is well known, minimum s, £-cut problems in undirected graphs can be 
represented by compact linear programs. Let G be a graph, let c : E(G) — > R+ 
be a capacity function, and let s and t two distinct nodes of G. For each node 
v G V^(G), introduce a variable z„, and for each edge e G E{G), a variable a; e . 
Then, the model reads: 



(16st) 

for all {u,v} EE(G). 

Hu [7] studied the following problem. Let G be a graph with "requirement" 
function c : E(G) — > H+ (say the data volume to be sent between two nodes of a 
network), and let K be the complete graph on V(G). A minimum-requirement 
tree is a spanning tree H of K minimizing 

Y, c(e)dist H (e), (17) 

eSB(G) 

where dist#(e) denotes the length of the path in H connecting the end nodes 
of e. 

For any edge / of H, define rjj(f) to be the capacity (=rcquirement) of the 
cut in G determined by the two components of H — f. This cut is called the 
fundamental cut induced by f. Then, (I17p is equal to 

E »•*(/)■ ( 18 ) 

feE(H) 

For a fixed spanning tree H of K and a fixed edge / of H, the capacity rjj(f) 
can be easily expressed as the optimal objective value of the linear program 



min Yl c(e)x e 




e£E(G) 




S.t. 


Zs = 0, 




Zt = 1, 


%{u,v} * %u ~ 


■*t, > 0, 


%{u,v} "T" %V 


z u > 



min ^ c {s)x( 

eS-E(G) 

s.t. 4 = 


for all v E U, 


4 = 1 


for all u G F(G) \ C/, 


4u,v } + 4 - 4 > o 


for all {u,w} G E(G), 
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where U is either of the components of H — f. This linear program can be 
expressed in terms of the compact spanning tree formulation (|10p by identifying 
z- and )JL- variables and fixing A- variables as follows: 



nin Yl c ( e ) x l 




eefi(G) 




•t. xl v} + M/,4,« - Hf,t,v > 0, 
X {u v} + M/,M ~ (if,t,u > 




for all {u,v}£E(G), 


A e = l 


for all e G J5(F), 


A e = 


for all e G £(A) \ £(#), 


(A,/z) satisfies (fTU)). 





where £ G / is fixed. Thus, fixing for each edge f oi H a. node £/ G /, (fT5|) is 
determined by 

min 5D 2 c ( e ) x l 

f£E(H)e£E(G) 
S.t. E^ „ } + M/,4/,« - Vf,tf,v > 0, 

s {u,i>} + /"/,4/,-u - /•*/,%,« > for all / e S(if), {u, v} G £(G), 

A e = 1 for all e G £(#), 

A e = for all e G A(A') \ £(#), 
(A, fj,) satisfies (fTU)l . 

Hence, it is obvious that the linear program 

mm Y^ Yl c ( e ) x i ( 19 ) 

feE(K)eeE(G) 
sX - X {u,v}+ Mtf^- M,t f ,v - ' 

x {u,v} + Vf,t f ,v - M/,t,,u > for all / G A(A), {u,u} G £(G), (20) 
(A, fj) satisfies (fTU)l . 

is a compact formulation for minimum-requirement trees. 

Theorem 3.1. Let G be a graph with capacity function c : E(G) — > R+, and 
let K be the complete graph on V(G). Moreover, let the polyhedron P be the 
set of all ({xf}f£E(K)> A, /i) satisfying (jTUJ) and (f2"0"l) . Then, for any vector 
(i x ' }feE(K):Kfi) G P, ({^''}/eB(ii'), A, /i) is a vertex of P minimizing (TIT))) 
i/ and cm/?/ i/ A is t/ie characteristic vector of a minimum-requirement tree H 
for (G,c), for each f G E{H), the vector (p>f,t f ,v)veV(K) * s the characteristic 
vector of the component of H — f containing tf, x' is the characteristic vector 
of the fundamental cut in G induced by f, and for each f G E(K) \ E(T), 
M/,«,u = x l — f or all u £ f, v € V. 

Proof. obviously is a lower bound of the objective values of the linear program 
to be considered, and hence, always an optimal vertex of P exists. 

Any vector ({x^ }f£E(K)> A, fi) G P is a vertex of P if and only if each 
component x, -, is chosen to be minimal, that is, 

x {u,v} = max {M/,t/,u ~~ M/,t/,«! M/,t/,« ~~ M/,t/,«} 
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and (A, /i) is a vertex of the extension II (constituted by (jTUJ)) of the spanning 
tree polytope PTree(-£0- So by Lemma [1.41 ({x'} f e E(K)i K 1^) 1S vertex if and 
only if satisfies the conditions mentioned in Theorem 13.11 

To conclude, for any vertex ({x '}/gB(jf), A,/i) of P and any spanning tree 
H of K, X = x E(H) implies E eGB ( G ) c ( e ) x l = r H(f) if / e E{H) and equals 
otherwise. Thus, ({x }/gB(if), A,/i) minimizes (j 19|) if and only if H is a 
minimum-requirement tree. □ 

Next, we turn to Gomory-Hu trees. Given a graph G with capacity function 
c : E{G) -> R + , a Gomory-Hu tree for (G, c) is a tree iJ with V(if) = V(G) such 
that for each edge e = {s,i} of H, the cut determined by the two components 
of H — e is a minimum s, i-cut of G. 

In the all-pairs minimum cut problem, one wants to find a minimum s,t-cut 
for all pairs s,t £ V(G). This problem can be represented in a compact form 
by Gomory-Hu trees. Let H be a Gomory-Hu tree for (G,c). For e = {s, t} £ 
E(K), define r({s, £}) as the capacity of a minimum s, i-cut of G. Then, 

r({w,w})= min r(e) for all {u, w} £ E(K), (21) 

e£E(P„„) 

where P u „ is the unique path in _ff connecting u and v. 

Hu [7] showed, based on an earlier result of Adolphson and Hu [I], that 
Gomory-Hu trees are minimum-requirement trees. Considering the proofs of 
the results in [TJ [7] it turns out that, conversely, minimum-requirement trees 
are Gomory-Hu trees. Below we restate Hu's result taking into consideration 
this equivalence. 

Theorem 3.2. Let G be a graph with capacity function c : E — > R + , and let K 
be the complete graph on V(G). Then, H C K is a Gomory-Hu tree for (G,c) 
if and only if H is a minimum-requirement tree for (G,c). 

Proof. Let H and H' be a Gomory-Hu and a minimum-requirement tree for 
(G,c), respectively. By definition of H, for each edge / = {s,t} G E(H) and 
each edge /' on the s, i-path in H' one has 

rw(f)>r H (f), (22) 

as rij(f) is the capacity of a minimum s,t-cut and the components of H' — f 
determine an s, t-cut. There exists a bijection ip : E(H) — > E(H') such that for 
each / = {s, t] £ E(H), <p(f) is an edge on the s, i-path in H', since H and H' 
are spanning trees. For details, see [1] or [T4J Section 15.4a]. So ((22]) implies 

E ^(/o= E r H'Mf))> E ^(/)- ( 23 ) 

feE(H') feE(H) feE(H) 

Since (fTTf and ([T%]) are equal, iJ' minimizes (|T^| , and hence all sums in (|2"3"|) are 
equal. This implies that i? is a minimum-requirement tree. Moreover, because 
of (|22), rj?'(/')) = r H ((y3 _1 (f )) for a U /' e £?(i?')- Hence, iJ' is a Gomory-Hu 
tree. D 

Corollary 3.3. Let G be a graph, and let c : E(G) — > R+ fee a capacity function. 
The linear program (|19p subject to (|10[) and (|20p is a compact formulation for 
Gomory-Hu trees for (G, c) . 
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The linear program (fT§j) subject to (flU)) and (f2"U|) does not exactly represent 
the two-level approach described in Section^ but it is in the same spirit. In this 
context, recall that (|21[) implies that H is a maximum spanning tree for (K,r). 
However, not every maximum spanning tree for (K, r) is a Gomory-Hu tree, see 
Schrijver [Mj Section 15.4]. Without this restriction, we only had to couple the 
min-cut formulations (|16stl) with a compact representation of the convex hull 
of the characteristic vectors of the edge complements of spanning trees (called 
co-trees), obtained, for instance, by a simple modification of system (I10[) . (In 
this context, note that Theorem 11.31 presumes to couple either maximization 
or minimization problems.) The restriction that not every maximum spanning 
tree for (K, r) is a Gomory-Hu tree leads then to the idea to identify [i- and 
z- variables as presented above. 

In what follows, we describe some consequences of Theorem 13.21 

3.1 The Gomory-Hu tree polytope 

Let G be a graph with capacity function c : E(G) — > R + . The Gomory-Hu 
tree polytope PGH-Tree(G, c) is the convex hull of the characteristic vectors of 
Gomory-Hu trees for (G,c). Using standard linear programming techniques, 
a compact formulation for PcH-Tree(G, c) can be derived from the linear pro- 
gram UnD subject to (fTni), (12"01) . 
Given a pair 

(P) min{w T x \ Ax > b} (JD) max{y T 6 | y T A — w T , y > 0} 

of dual linear programs. If both (P) and (D) have feasible solutions, then 

Q := {{x, y) | Ax > b, y T A = w T , w T x -y T b = 0,y> 0} 

is a polyhedron that consists of all vectors (x, y) such that x and y are optimal 
for (P) and (D), respectively. 

If (P) is the linear program (TK)|) subject to (ITU]) . (|2T)|) . then we conclude from 
Theorem 13.21 that Pgh-Tioc(G, c) is the projection of Q onto the A-space. 

Corollary 3.4. Let G be a graph, and let c : E(G) — > R+ be a capacity func- 
tion. Then, there exists an 0{\V(G)\ A + \V(G)\ 2 (c)) extended formulation of 
PGH-TVee(G, c) , where (c) denotes the input length of c. 

3.2 The minimum T-cut problem 

Let G be a graph with capacities on the edges c : E(G) — > R+, and let T be 
a subset of the nodes of G of even size. A T-cut is a cut 5(U) of G such that 
\T n U\ is odd. In the minimum T-cut problem one wants to find a T-cut S(U) 
of G minimizing c(8(U)). The polyhedral counterpart to the minimum T-cut 
problem is the T-cut polyhedron, defined as the set of all x £ R- 6 ^ ) for which 
there exists a convex combination y of the characteristic vectors of T-cuts such 
that x > y. 

The minimum T-cut problem can be solved in polynomial time. Various 
algorithms are available and, among them, the famous algorithm of Padbcrg 
and Rao [12] that computes a Gomory-Hu tree for (G, c) , and selects among 
the fundamental cuts a T-cut of minimum capacity. This T-cut has minimum 
capacity among all T-cuts of G. 
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This algorithm can be used as as basis to derive a compact formulation for 
the minimum T-cut problem. However, the inequalities of our formulation will 
depend on the input c, and as a consequence, this is only a formulation for the 
minimum T-cut problem but not for the T-cut polyhedron. 

Let K be the complete graph on V(G), let H be a Gomory-Hu tree for (G, c), 
and let A its characteristic vector. Recall that for / G E(H), the capacity of 
the fundamental cut induced by / is the optimal objective value of the linear 
program 



min 2_, c i e ) x e 






(24f) 


e£E(G) 








s.t. A e = 1 




for all e G E(H), 


(25) 


A e = 




for all e e E(K)\E(H), 


(26) 


X {u,v} + Vf,tf,u - 
f 


- tJ-f,t f ,v > 0, 

- Hf,t f ,u > 


for all {u,v} eE(G), 


(27f) 



(A,/z) satisfies (flOl) . 

where £/ G / is fixed, while for / G E{K) \ E{H) the optimal value is zero. 
Introduce for each / G E(K) a (binary) variable Vf such that Vf = 1 if and only 
if / G T(T/) and the fundamental cut induced by / is a T-cut. Denote the edge 
set corresponding to T-cuts by F. Then, the following mathematical program 
is a disjunctive programming approach for finding a minimum T-cut among the 
fundamental cuts of H: 

min E E < e )y f e ( 24 ) 

feE(K)e<EE(G) 

s.t. (A, n) satisfies (ITU | . ([2"5]) . (|2"o) . 

({i'}/6E(Jf), AO satisfies ((273) for all / e E(K), (27) 

/ 1 if/eT 

•* [0 otherwise, 

^ ^ = 1, < $ < v, (29) 

feE(K) 
y{ >xt + f -l, y{ > for all e e T(G), / £ E{K). (30) 

To see the correctness of this formulation, let {{x* tD^} f£E(K)i A, /x, ;/, i?) be 
a vertex of the polyhedron determined by the inequalities of this formulation. 
Recall that A is the characteristic vector of H. Moreover, if / = {s,t} G E(H), 
fif yS ^,fif y t,* are the characteristic vectors of the components of H — f, and x' 
represents the associated fundamental cut. Otherwise, that is, in case A/ = 0, 
X? = /J.f, s .* = M/,t,* = 0- Moreover, since ({x-^, y^}f£E(K)> A, /i, v, -&) is a vertex, 
inequalities (f2"5|) . (f2U|) imply that d is a unit vector with -d g = 1 for some g <E F. 
For the same reason, ({x-' 7 , yf}f<£E{K), A, ^, ^, ^) satisfies at least one of the two 
equations in (|5P|) at equality for each pair of edges e G E(G),f G E(K). Hence, 
y* = for all / G T(if) \ {5}, while y 9 is the characteristic vector of a T-cut. 

In an instance of the minimum T-cut problem neither H nor F are part of 
the input. Therefore, to derive a compact linear formulation from the above 
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mathematical program we have to remove the fixing equations (J2"H) , ([25)1 and to 
replace the fixing equations ([25)1 by an appropriate system of inequalities. 

In a first step, let us assume that H but only H is part of the input. One 
way to determine v is as follows. Fix any node r G T, and for each s£T\ {r}, 
denote by P s the unique r, s-path in H. Then, for any edge / of H, f G F if 
and only if the number of paths P s using / is odd. This, in turn, is the case if 
and only if the symmetric difference of the paths P s contain /. 

Let s G T\ {r}. Introducing for each edge e G E{K) a variable 7r|, one easily 
checks that the characteristic vector of P s is, for instance, determined by the 
system 

y°(8(r)) = y s (6(s)) = 1, 

y s (S(v) \ e) - y% > for all v £ V(iQ, e G <J(v), (31s) 

< y s < A. 

Next, we have to express the symmetric difference of the paths P s by linear 
inequalities. As it is well-known, the characteristic vector of the symmetric 
difference XAY := (lUf)\(Ifl Y) of any sets X,Y C E(K) is determined 
by the system 



"/ < Xf + X Y f 




for all / G E{K), 


u f >xj- x] 




for all / G E(K), 


a f > X Y + xf 




for all / G E(K), 


a f <2- X f- 


xj 


for all f EE(K). 



Since the A-operator is associative, the symmetric difference of the paths P a 
can be determined sequentially. For this, let (so, S\, . . . , Sk) be any order of the 
nodes in T\{r}. Then, defining Dq := P SQ and Di := Di_iAP Si for i = 1, . . . , k, 
it follows that F = Dk- Thus, introducing for each step i G {0, 1, . . . , fc}, an edge 
variable vector a 1 to represent Di, we see that equations (|!?5j) can be replaced 
by the system 



y Si , A) satisfies 


dsg 


a - y S0 


= 


a 1 + a l ~ l +y Si 


>o 


a 1 -a l ~ l +y Sl 


>o 


a 1 + a l ~ l -y Si 


>o 


a 1 + a l ~ l +y Sl 


< 2 


a k - v 


= 



, i = 0,l,.. 


• j ^j 


(31) 

(32) 


for i = 1 , . . 


■ ? ">, 


(33) 


for i = 1 , . . 


■ ? «) 


(34) 


for i = 1 , . . 


■ ? "-) 


(35) 


for z = 1 , . . 


■ ? "<) 


(36) 
(37) 



In the remainder of this section we consider the linear program (1241) subject 
to Uni), d25j)-(EZI), (EH)-©- Removing the fixing variables ([25]), &§§, we obtain, 
of course, a compact formulation for a relaxation of the minimum T-cut problem. 
This implies that the system (fTU|> . (|2T|) . ([2^1) - ([3T)l characterizes a polyhedron 
that contains the T-cut polyhedron, but we conjecture that this system does 
not determine the T-cut polyhedron. However, using standard techniques, a 
compact formulation for a particular instance of the minimum T-cut problem 
can be derived as follows. By Corollary 13.41 there exists a compact extension 
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Q of the Gomory-Hu polytope P GH-Tree{G , c) . We assume that Q is given in 
A, [i, if- space. Let (/i, / 2 , . . . , f m ) be any ordering of the edges of K. Then, the 
unique optimal solution of the linear program 

m 

min ^2 J - 1 A ei 

.7=1 
s.t. (A, /i, 77) 6 Q 

is the characteristic vector of the minimal lexicographic Gomory-Hu tree with 
respect to the ordering. Q is determined by a compact system of the form 

A >b, 

and by the construction in Subsection 13. 1[ we may assume that this system con- 
tains (ITU|) as subsystem. The characteristic vector of the minimal lexicographic 
Gomory-Hu tree is the projection of the polytope determined by 

A I At >b, ir T A = d T , d T J At J - TT T b = 0, tt > (38) 




onto the A-space, where d T := (1, 2,4, . . . , 2 m_1 ,0, 0, . . . , 0), and hence, the 
linear program ([24]) subject to (l27l) . (l29l) - ([38|) is a compact formulation for the 
minimum T-cut problem. However, we note that (I38p contains inequalities 
some of whose coefficients have input length |2?(lif)|, and more important, some 
depend on c, and thus this is not a formulation for the T-cut polyhedron. 

4 Conclusion 

Motivated by a well-known result of Balas on disjunctive programming, we stud- 
ied extended formulations in connection with a simple two-level optimization 
scheme that also can be derived as a special case from the more general frame- 
work branched polyhedral systems of Kaibel and Loos [5]. Using this scheme, we 
gave compact formulations for the spanning tree approximation of the Stcincr 
tree problem and for Gomory-Hu trees. Using the Gomory-Hu tree formulation, 
we also derived a compact formulation for the minimum T-cut problem whose 
inequalities, however, depend on the objective function. 

A very interesting question in this context is if the compact formulation for 
the minimum T-cut problem can be modified in such a way that it extends 
to a compact formulation for the T-cut polyhedron. However, answering this 
question seems to be a complicated undertaking, since the separation problem 
for the perfect matching polytope defined on a graph G can be modeled as linear 
optimization problem over the V(G)-cut polyhedron. Such a formulation would 
also imply a compact formulation for the perfect matching polytope - a long 
standing unsolved task. 

A powerful theory on coupling extended formulations has been developed, 
especially due to the contribution of Kaibel and Loos [5] . In future research, we 
are interested to find further examples, where the coupling of extended formu- 
lations yields to compact formulations for polynomially solvable combinatorial 
optimization problems. 
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